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Abstract 


I.  A  fornmla  is  derived  which  expresses  the  perturbed  scattering 
amplitudes  of  a  combination  of  two  arbitrary  cylinders  as  a 
function  of  the  iinperturbed  scattering  amplitudes  of  the  in- 
dividual cylinders.  The  formula  is  valid  when  the  spacing  of 
the  scatterers  is  large  compared  to  their  dimensions.  It  in- 
volves derivatives  of  the  scattering  amplitudes  with  respect 
to  the  angles  of  incidence  and  of  ohyersation.  Interaction 
terms  of  degrees  d"-'-/^,d"  ,  and  d~  '   are  taken  into  account, 
where  d  is  the  spacing.   Verification  is  obtained  in  a  special 
case.  The  result  is  employed  to  calculate  the  total  scattering 
cross  section. 

II.  The  method  of  I  is  extended  to  cover  the  three  dimensional  scalar 
problem  for  two  bodies  of  arbitrary  shape.  All  interaction  terms 
of  order  d   and  d   are  given. 


Table  of  Contents 

Page 

I.  Two  Dimensional  Case  1 

1.  Introduction  1 

2.  A  statement  of  the  problem,  and  of  the  method  of 
analysis  3 

3.  Expansion  of  the  scattered  wave  in  terras  of  plane 

waves  6 

A*  Expansion  of  the  response  of  cylinder  A  in  a 

neighborhood  of  cylinder  B  6 

B.  Expansion  in  a  neighborhood  of  A  for  the  wave 

scattered  by  B  10 

C»  Expression  of  the  scattered  waves  in  terms  of 

plane  waves  10 

D,  Elimination  of  f,  12 

h»   Calculation  of  the   interaction  13 

5.  A  special  case.  Scattering  by  two  conduction 
circular  cylinders  18 

6.  An  explanation  of  the  result  2U 

7.  The  total  scattering  cross  section  of  a  combination 

of  two  identical  circular  cylinders  26 

8.  Additional  remarks  28 

II.  Three  Dimensional  Case  30 
1.  Statement  of  the  problem  30 
2«  Expansion  of  the  scattered  ;rave  in  terms  of  plane  waves33 
3.  Calculation  of  the  perturbed  far  field  amplitude  37 
h»  Explanation  of  the  result  39 
References  ill 


-  1  - 

I.  Two  Dimensional  Case 
1.  Introduction 

The  present  paper  deals  with  the  diffraction  of  plane  electromagnetic 

or  acoustic  waves  by  a  pair  of  parallel  cylinders  of  arbitrary  shape. 

The  diffraction  by  each  cylinder,  in  isolation,  is  assumed  known,  and  the 

diffraction  by  the  configuration  is  calculated  explicitly,  in  terms 

of  these  data.  An  approximation  is  involved  which  will  be  described 

below. 

The  question  of  multiple  scattering  has  already  been  treated,  but 
in  less  detail,  by  a  number  of  writers.  A  brief  sketch  of  the  history  of 
the  problem  follows,  with  emphasis  on  those  treatments  of  the  problem 
whose  accuracy  increases  with  the  spacing.  We  recall  Reiche  and  Schafer"-  -' 
who  were  the  first  to  have  given  a  wave  theoretical  discussion  of  the 
finite  grating  of  circular  cylinders.  These  authors  neglected  the 
interaction  between  the  cylinders}  their  work  was  therefore  valid,  in 
principle,  in  the  limit  of  large  spacing.  A  very  general  expression  for 
the  diffraction  by  an  arbitrary  assemblage  of  circular  cylinders  was 
given  by  Twersky  ^'^,  who  took  all  orders  of  interaction  into  account. 
The  method  depended  heavily  upon  the  separability  of  the  circular 
geometry,  and  the  most  general  form  of  the  result  was  too  complicated  to 
be  discussed.  However,  Twersky  found  that,  if  he  proceeded  to  the  limit 
of  large  spacing,  he  could  simplify  his  res\ilt  immensely  and  achieve  a 

perspic\x)us  discussion  of  the  correction  to  single  scattering.  This  woi^c 

-1/2     -1 
gave  a  correct  accoimt  of  terms  of  degree  d    and  d   in  the  spacing. 

This  success  led  to  a  general  investigation  of  the  large  spacing 
approximation  by  Karp  '^-',  who  showed  that  for  cylinders  of  arbitrary 
shape,  the  leading  terms  of  the  interaction  correction  could  themselves 
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be  expressed  explicitly  in  terms  of  non-interaction  or  single  scattering 
results,  the  latter  being  regarded  as  given.  In  fact  the  interaction  term 
could  be  regarded  as  being  composed  of  the  response  of  each  cylinder  to 
a  plane  wave  arriving  from  the  direction  of  the  other  cylinders.  The 
techniques  of  [li]  were  exploited  by  Karp  and  Radlow  in  the  analysis  of 
a  grating  of  cylinders.  Similar  methods  were  used  by  Karp  and  Russek^  -' 
in  expressing  the  approximate  solution  to  the  problem  of  diffraction  by 
a  wide  slit  in  terms  of  the  well  known  solution  for  the  half  plane  problem. 

The  purpose  of  the  present  paper  is  to  extend  the  work  of  Karp  ^-' 
so  as  to  take  into  account  higher  order  terms.  JT:ist  as  in  [3_| ,  the 
cylinders  are  arbitrary,  arxi  the  scattering  by  each  cylinder  in  isolation 
is  assumed  as  given.  But,  it  is  found   that  even  the  higher  order 
correction  terms  can  be  calculated  generally,  simply  and  explicitly  in 
terms  of  the  single  scattering  data  used  in  [3]  for  calculation  of  the 
leading  terms.  This  is  the  principal  result  of  the  present  paper. 

The  general  calculation  was  carried  out  so  as  to  include  all  effects 

-I/2   -1     -'?/2 
of  order  of  magnitude  d    ,  d   and  d  ^'  ,  where  d  is  the  spacing  of 

the  cylinders.  For  purposes  of  comparison,  Twersky's  calculation  for  a 

pair  of  circular  cylinders  was  continued  so  as  to  include  terms  of  order 

d  '  }  this  special  calculation  by  the  repeated  application  of  addition 

theorems  was  then  shown  to  agree  with  the  general  result  of  the  present  work, 

when  the  latter  result  is  specialized  to  the  case  of  circular  cylinders. 

The  result  is  used  to  calculate  the  total  cross  section  for  a  pair  of 

circular  cylinders  as  a  function  of  the  spacing  and  the  known  unperturbed 

(or  non-interaction)  scattering  amplitude  functions. 

See  also  [5] 
**Cf.  [7] 
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2,     A  Statement  of  the  Probleirij  and  of  the  Method  of  Analysis. 

We  would  like  to  know  the  scattering  pattern  of  a  combination  of 
two  infinite  parallel  cylinders  in  terras  of  the  scattering  patterns 
which  these  cylinders  would  have  if  they  were  isolated  from  each  other. 
In  other  words,  we  want  to  obtain  a  functional  relationship  between  the 
imperturbed  and  the  perturbed  scattering  patterns  of  the  cylinders.     Such 
a  relation  is  desirable  because  it  sijiiplif  ies  the  calculation  of  the 
pattern  for  the  combination.     If  we  can  calculate  the  unperturbed  patterns, 
we  need  only  insert  them  in  this  relation  to  obtain  the  perturbed  patterns. 
The  relation  is  useful,  moreover,  even  if  the  shapes  of  the  cylinders 
are  so  complicated  that  we  cannot  separate  variables  or  if  calculation 
by  separation  of  variables  is  too  tedious   .     Also  in  such  cases,   the 
unperturbed  patterns  might  be  measured  experimentally  for  all  angles  of 
observation  and  these  unperturbed  patterns  might  then  be  substituted 
into  the  relation  obtained  here  to  yield  the  perturbed  patterns. 

The  situation  is  the  following:     A  plane  wave  of  unit  amplitude  is 
incident  upon  the  two  parallel  cylinders  A  and  B  (Figure  1). 

LEFT  CYLINDER  A  RIGHT  CYLINDER  B 


Incident 

Plane 

Wave 

FIGURE  1 


See   [k]  for  the  use  of  an  even  less  accurate  relation  for  theoretical 
purposes . 
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To  avoid  ambiguity  in  the  definition  of  the  spacing  d,  we  define 

a  coordinate  system  for  each  cylinder.  Let  A*  and  B  be  circular 

cylinders  circumscribed  about  A  and  B  respectively.  Let  a  and  b  be 

the  respective  radii  of  A*  and  B*.  We  shall  let  Z  be  the  axis  of  A* 

a 

and  Z.    be  the  axis  of  B   .     The  problem  is  two-dimensional,  and  we  shall 
operate  in  a  plane  perpendicular  to  the  Z  axes.     The  respective  intersections 
of  this  plane  with  the  Z     and  Z,    axes  will  then  be  the  origins  of  the 
coordinate  systems  of  A  and  B   .      We  can  now  define  the  spacing  d  as  the 
distance  between  the  axes  of  A     and  b'-  that  is,  the  distance  between 
the  two  origins  (Figure  2). 


FIGURE  2 

We  make  the  following  assumptions: 

1)  The  individual  complex  scattering  pattern  is  known  when  each  cylinder 
stands  alone  in  space. 

2)  d  >  >  X  where  X  is  the  wave  length  of  the  incident  wave.* 

3)  d  >  >  a  and  d  >  >  b. 


Numerical  calculation  in   [$] ,  which  are  based  on  a  less  accurate  procedure, 
showed  that  d  could  be  as  small  as  one  wavelength  without  impairing 
the  accuracy  materially.     The  present  result  would  therefore  allow  an 
even  smaller  spacing,  provided  the  cylinders  are  kept  sufficiently  small 
compared  to  the  spacing. 
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Our  object  is  to  find  a  functional  relation  between  the   scattering 
pattern  of  the  combination  and  those  of  the  isolated  component  cylinders. 
This  cannot  be  accomplished  by  simple  superposition  of  the  unperturbed 
patterns  of  the  components  since  the  individual  pattern  of  each 
component  cylinder  is  modified  by  the  field  scattered  by  the  other 
cylinder.    We  must,  therefore,  consider  the  interaction. 

We  shall  asstime  that  the  response  U    of  each  cylinder  to  a  plane 

s 

wave  is  of  the  fonn 


ikr  00  f  (9  &  ) 

(1)  u    =  -^ —  Y~  -2 2_  for  large  r 

S    ,r—  ^—        n 

/r  n=o      r 


where  r  is  the  distance  from  the  axis  of  the  circumscribed  circular  cylinder, 

9  is  the  angle  of  observation,  and  &  is  the  angle  of  incidence.  Both  Q 

and  ©  are  measured  from  the  x  axis  of  each  cylinder.  The  x  axes  are 
o 

collinear  with  d. 

As  will  be  explained  below,  at  a  sufficiently  large  distance  from 

the  cylinder,  the  scattered  field  U  resembles  a  plane  wave.  This 

s 

approximate  plane  wave  elicits  a  response  from  the  other  cylinder, 
perturbing  its  scattered  field.  This  response  also  has  the  form  (1)  and 
in  turn  perturbs  the  field  scattered  by  the  first  cylinder.  We  can  carry 
out  successive  calculations  for  this  process  until  the  desired  order  of 
accuracy  is  obtained.  When  the  perturbed  patterns  have  been  calculated, 
they  can  be  superposed. 

We  shall  deal,  in  this  paper,  with  interaction  terms  of  degrees 

-1/2   -1     -3/2 
d    ,  d   and  d    .  The  procedure  involves  a  new  kind  of  expansion  of 
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the  waves  scattered  by  each  cylinder  about  the  origin  located  in  the 

-3/2 
other  cylinder.     In  order  to  insure  that  all  terras  up  to  order  d 

are  contained  in  the  result,  we  niust  include  them  in  the  first  expansion. 

We  then  find  that  the  field  scattered  by  a  given  cylinder  can  be 

represented,  in  the  neighborhood  of  the  other  cylinder,  aa  a  plane  wave, 

plus  additional  terms  which  are  recognized  as  derivatives  of  a  plane 

wave  with  respect  to  its  angle  of  incidence.     The  simple  way  of 

expressing  the  higher  order  excitiations  is  what  enables  ua  to  calculate 

the  higher  order  responses  conveniently, 

3.     Expansion  of  the  scattered  waves  in  terms  of  plane  waves* 

A,     Expansion  of  the  response  of  cylinder  A  in  a  neighborhood  of 
cylinder  B. 


Let  us  consider  what  hiappens  when  a  plane  wave  is  incident  upon  A, 
The  wave  function  for  a  plane  wave  is  (Figure  3). 

ik(x,co3  ©     +  y    sin  6  ) 
(2)  U^  -  e        *  °       ^ 

The  wave  scattered  by  A  in  response  to  the  plane  wave  is  represented 
by  an  asjTnptotic  solution  of  the  reduced  wave  eqviation. 
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(3)  D.E. 


(^  +  k^)  U  =  0  . 


The  boundary  conditions  are  such  that  the  diffraction  problem  for 
each  cylinder,  and  for  the  combination,  is  well  posed.  They  may  be, 
for  example, 


(U)     B.C. 


au 


a)       U  «  0 
or    b)    |2.0 


I 


dU 


or    c)     CU  +  D|ii  -  0 


where  x-  is  the  nomjal  derivative  of  U  and  C  and  D  are  constants. 
dn 

Alternatively,  one  or  both  of  the  cylinders  may  be  filled  with  dielectric 
materials. 

The  radiation  condition  for  an  outgoing  wave  i.e.; 


(5) 


lim 
r->oo     r 


^^Ij 


eu 

8r 


-  ikU    -  0 


is  imposed  on  all  scattered  fields  which  occur. 
The  solution  of  (3)  has  the  form 


(6) 


u  -  u  +  u 

i    a 


where  U.  is  the  incident  field  and  U,  is  the  field  scattered  by  cylinder  A. 
We  assume  that  U  may  be  represented  in  the  asymptotic  form 


(7) 


U.  - 


ikr 


'  a 


f  °  1 

'(9  ,9  )+  -L- (9  ,9  )+  0(r  "^) 
a*  o    r    a'  o    ^  a   I 

a  J 


The  letter  "  a"  signifies  that  the  variable  in  qiestion  refers  to  cylinder 
A,  and  the  superscript  "  -"  signifies  that  the  pattern  is  unperturbed. 
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If,  as  we  have  aasvuned,  B  is  sufficiently  small  in  relation  to  its 
distance  from  A,  then  the  wave  scattered  by  A  is  practically  a  plane 
wave  in  the  neighborhood  of  B,  and  we  may  imagine  such  a  wave  incident 
on  B,  We  can  demonstrate  an  explicit  representation  of  this  approximately 
plane  wave  by  expanding  U  in  a  neighborhood  of  B.  We  shall  express 

Si 


the  expansion  in  powers  of  d 


-1/2 


in  a  rectangular  coordinate  system 


with  origin  at  the  center  of  B  and  carry  out  the  calculation  up  to 


order  d 
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The  calculation  proceeds  as  follows: 


^l<-^ ^ 


FIGURE  U 


Let  P  be  a  point  in  the  neighborhood  of  B. (Figure  U).     Let  r  =  the  distance 
from  the  axis  of  A*   to  P.     We  see  from  Figure  U  that  r     =  /(d+x^^)  +  y,    and 
we  shall  henceforth  omit  the  subscripts   "  a"    and  "  b"  as  a  matter  of 
convenience.     Then 

(8)     e""    -    e^<'^*'='[l  *  ^  *     "^'-  ^A  "V*    ♦    ff(d-3)l 


(«    r-^^  .  Cd=.  2^  .  .^  .  .^)-^A  .k^.-^ 


^     ♦tfCd-^/^) 


(M)  r-3/2  .  d-3/2[l  -  3  I  *  ©-(d-^)]  .  -1^  *  <y(d-5/2) 
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(IDa)    ®  -  ®o  °  arctan  ^fe  "  d  "  ^  *  ^"^^"^^ 


(10b)  fn°(»,%)  -  f 


r(°'^o^^[vr(°>%)}4'^'^°^''^°^] 


where  D^,  2  ^  .  Using  (liDa),  we  can  rewrite  (lOb)  in  the  form 


2d 


Substituting  (8)  -  (11)   iato   (?),  we  obtain 


(12) 


1*^*  W-  lA  i^V  .a(d-3) 


t 


1  ^      +fVrH-^/2> 


> 


^.  -  e-'-'[^ 


im 


We  have  thus  obtained  the  field  scattered  by  A  as  it  appears  in  a 
neighborhood  of  B.     We  must  now  carry  out  the  corresponding  expansion  for 
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the  field  initially  scattered  by  B,  Note  that  the  leading  term  of  (12) 
is  a  plane  wave,  i.e.,  a  constant  multiple  of  e    ^  as  we  explained  earlier. 
B.  Expansion  in  a  neighborhood  of  A  for  the  wave  scattered  by  B« 

The  same  type  of  procedure  as  that  used  above  gives  us  the 
expansion  in  a  neighborhDod  of  A  of  the  wave  initially  scattered  by 
B.  The  expression  here  corresponding  to  (12)  is: 

(.3)   ^,.e-(-^)[f>;^ 

l/2(x  *  ikyh   f^("*»o^  -  y^e^r^"'^^  *  ^r("'^) 


7^ 


*  e'(d-^/2) 


] 


The  unperturbed  fields  scattered  by  either  cylinder,  (12)  or  (13) » 
are  the  excitations  of  the  other  cylinder.  We  assumed,  to  begin  with, 
that  we  knew  the  unperturbed  response  of  each  cylinder  to  a  plane  wave 

excitation  for  all  angles  of  incidence.  If  (12)  and  (13)  were  plane 

■it- 
waves,  we  could  calculate  the  responses  for  the  second  scattering, 

(12)  and  (13),  however,  are  not  plane  waves.  But  this  impediitient  does 

not  prevent  the  calculation  of  the  effect  of  further  scattering.  The 

reason  is  that  (12)  and  (13)  may  be  represented  in  terms  of  plane  waves 

by  appropriate  substitutions.  This  representation, which  will  enable  us 

to  calculate  successive  scattering,  will  now  be  given. 

C.  Expression  of  the  scattered  waves  in  terms  of  plane  waves. 

We  can  reduce  the  further  scattering  of  singly  scattered  waves  to 
the  scattering  of  plane  waves  by  expressing  (12)  and  (13)  in  terms  of 


This  WDuld  be  the  method  used  in  [3]  . 
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plane  waves  and  derivatives  of  plane  waves.  Noting  that  a  plane  wave 

is  represented  by 

ik(x  cos  &„  ■*•  y  sin  9  ) 
(Hi)         v(e^)  =  e         ° 

we  observe  that 


(15)  ikye^^  -    v^  (O) 

o 


(16)       -ikye~  "    "''^B  ^"^ 

o 


(17)  ik(-x*-iicy^)e^^  =     ^e  9  (O) 

o  0 


(I8)ik(x+iky2)e"^^     =     v^  ^  (n) 

o  o 


Substitution  of  (1$)  and  (1?)  in  (12)  and  of  (16)  and  (18)  in  (13) 
yields  the  following  representation  of  the  scattered  fields  in  terms  of 
plane  waves: 

•,,^rv(o)f^°(0,9   ) 
/-qv  .x>         ikd       ^    ^  o   '    »  o 

^^9)  K  =  «      [ ^T^ 

^[d|^v(o)]  ff  (0,9^)4^  [d^^v(o)]  D^f-(o,e^)  .  v(o)  f-(o,e^: 


d5/2 


..^,     v(T,)f^°(TT,&     ) 

roni  Ti°   -  ^^d               o        *  0 
(20)  u;  -  e  375 


^K"^"^]  'o°("^^)*  ^[v^")]  Vo°("*^o^^(")  ^i°("'^: 


) 
. ^^ ^d-^/2) 


7372 
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D»  Elimination  of  f-, 

We  note  that  the  numerators  of  (19)  and  (20)  are  sums  of  terms, 
consisting  of  plane  waves  and  their  derivatives,  namely  the  v's,  and 
coefficients  which  are  independent  of  9,  namely  the  f 's.  We  see  that 
these  formulas  are  expressed  in  terms  of  both  f  and  f,.  (For  the 
meaning  of  f  and  f,  see  (?)  above).  An  advantage  would  result  from 
the  elimination  of  f, ,  since  we  could  then  express  the  result  in  terms 
of  the  scattering  amplitude  of  the  far  field  without  having  to  know 
the  scattering  amplitudes  of  further  asymptotic  terms. 

We  eliminate  f,  by  expressing  it  in  terms  of  f  .  This  can  be 
done  by  means  of  a  recursion  formula.  The  recursion  is  obtained  by 
substitution  intn  (3)  of  the  assumed  asymptotic  representation  (?), 

of  pny  radiating  solution  o-f  the  reduced  wave  equation.  When  we  equate 
the  corresponding  inverse  powers  of  r,  we  find  that 

This  recursion  is  useful,  also,  for  calculations  of  higher  degree 

than  we  are  considering  here.  Since  we  want  to  express  f,  in  terms  of 

f  ,  we  need  use  it  only  for  the  value  n  =  1. 
o 


(22)         h'm 


■P  f  +  D^f 
U  o    9  o 


This  is  a  different  type  of  recursion  from  that  obtained  for  large  k  by 
Keller,  Lewis,   and  Seckler'--'   although  it  is  similar  in  form. 


The  representation  is  asymptotic  for  large  r,  9  and  k  being  held  fxxed. 
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If  we  now  substitute  (22)  into  (19)  and  (20),  we  obtain; 


(23)     U°  =  e 


ikd 


v(o)f;"(o,e^) 

— :;m — 


1 

2ik 


d;  v(0) 
L  o 


L  o    . 


Vo"(°'^o^^4-^^) 


^f^°(o,e  )+D^f^°(o,e  ) 

U  o   '  0   a  0   '  o 


,3/2 


.^(d-5/2) 


(2U)     U^  =  e 


ikd 


v(0)f^°(n,e^) 

— ::i/2 — 


1 
231 


d|  v(n) 
L  0 


f^°( 
o 


»'»o'*i[v<")]v?("'»o)*  SS'*"'!  -'^('".*o>*'>|f^('"«o>] 


•d(d-5/^) 


The  above  expressions  (23)  arkd  (2U)  represent  the  responses  of  the 
cylinders  to  the  original  incident  plane  wave  and  these  responses  are  given 
near  the  other  cylinder  in  terms  of  plane  waves.   (23)  and  (2U)  are  also 
excitations  for  the  second  scattering,  Succesaive  application  of  these 
formulas  will  yield  the  desired  degree  of  interaction. 

U.  Calculation  of  the  Interaction. 

We  have  expressed  in  (23)  and  (2U)  the  fields  singly  scattered  by 
A  and  by  B,  as  they  appear  in  a  neighborhood  of  the  second  scatterer  and 
have  moreover  expressed  them  in  terms  of  plane  waves  and  derivatives  of 


-  Hi  - 

plane  waves.  We  may  now  imagine  this  combination  of  plane  waves  and 
their  derivatives  to  be  incident  upon  the  second  scatterer.  The  linearity 
of  these  expressions  enables  us  to  say  that  the  responses  of  A  and  B  to 
incident  derivatives  of  plane  waves  are  equal  to  the  derivatives  of  the 
responses  of  A  and  B  to  the  incident  plane  waves.  Since  we  already 
know,  by  assumption,  the  unperturbed  responses  of  A   and  B  to  an  incident 
plane  wave,  we  have  reduced  the  second  scattering  to  the  previous  case, 
namely  the  first  scattering. 

If  we  carry  out  this  process  to  the  extent  of  three  successive 

-1/2   -1 
scatterings,  we  can  obtain  interaction  terms  of  degrees  d    ,  d   and 

-3/2  -I/2 

d  '  ,  The  terms  of  degree  d     result  from  double  scattering, 

those  of  degree  d~  result  from  triple  scattering,  and  those  of  degree 

-3/2 

d     result  partly  from  double  and  partly  from  quadruple  scattering. 

For  the  sake  of  simplicity,  the  following  results  will  be  expressed 
in  a  normalized  coordinate  system  (Figure  5)  with  a  common  origin  midway 
between  the  origins  of  the  coordinate  systems  located  in  the  scatterers. 
As  a  matter  of  convenience,  we  shall  omit  the  subscript  zero  from  the  f 's. 

Q 


FIGURE  5 
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Let  a  be  the  angle  of  incidence  of  the  original  plane  wave.     The 

perturbed  patterns  will  then  be 

.,d 
-ik— cos  a 

(25)  f^(e,a)   =  e       2  f^°(e,a) 

ikrrcos  a   .,  ,  , 
,   e     ^  e^^f^(n,a)f^°(e3n) 

-ik?rco3  ct   .,  „  ,  , 

^  e       ^  e^^V°(0,a)f^°(n,0)f^°(e,Ti) 

d 

^'°'  °  ik3d    ^ 
+  5 5^^-^ f^°(TT,a)f^°(o,n)f^°(n,o)f^°(»,n) 

ik|co3  a  ii^^/ 
*-2t75 ^—  f^°(.,a)D2  f^°(0,Tt)  .  2D.f^^n,a)D.  f^°(»,n) 


T^f^°(n,a)  .  D|f^(n,G)|  f^°(©,n) 


^e)(d-5/2) 


and  ,  .,  d 

i4cos  a  ^                  -^?°°^  "  ikd 
(26)  f^(9,a)  =  e    ^         f^(»,a)  +5 j^-^—     f^°(d,a)f^(9,o) 

^1^°^  ^  ik2d  ^ 

+  ^ J— 2 f^(Ti,a)f^(o,n)f^°(e,o) 

-ik|co3  a 

e  e  -ao/        \o,bO/        \«ao/        \-.bo. 


*  ' j^ f**"(o,c)f'-(n,o)f^"(o,n)f""(9,0) 

-ikfcos  a  ^  / 

*  ^       ^3/^     "         f  f^°(o,a)D^^f^O,o).  2Dgf^°(o,a)D^  f^°(S,o) 

+  |'^f^°(o,a)  .  D|f^°(o,a)l  f^(ft,o: 
-©'(d-^/2) 
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We  nate  that  the  successive  powers  of  d  in  (25)  and  (26) 

represent  the  various  degrees  of  interaction.     The  scattered  far  fields 

resulting  from  the  interaction  will  have  the  form 

ik(r  +  I  cos  9) 

(27)  U     =  ^ f^(e,a) 

I/? 

and 

ik(r  -  I  cos  e) 

(28)  U     ».2 f^(&,a) 

|/r 

where  f^(^,a)  and  f^(&,a)  are  the  perturbed  scattering  amplitudes 
given  by  (2$)   and  (26).     The  sum  of  the  scattered  fields  has  the  following 
far  field  representation  which  is  the  sum  of  (2?)  and  (28). 


Ucr     nr-       -if 
(29)  U     "  U^  +  U^  »  ^  /  ^     e     ^  F(»,a) 


where  F(»,a)  is  the  scattering  amplitude  of  the  combination  of 
cylinders.     We  see  by  comparison  with  (2?)  and  (28)   that 


,-T.     i^I  /  i4  cos  e  -ik|cos  e 

(30)     F(#,a)-/^e^    e     ^  f^(^,a)     .     e       ^  ^i^^^^) 
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A  nt>re  explicit  formula  for  the  far  field  is  obtained  by  combining   (25), 
(26),    (29)   and   (30).      The  far  field  can  then  be  written  as  follows: 

ikr        il<7T(  cos  &  -  cos  a )  a 


^  e 


(31)  u  S  ^__  .  e 


l^ 


-ikT?(cos  9-  cos  ot) 
f'°(&,a)  *  e       2  f^°(e,c) 


ikd 


ik'K(cos  &  +  cos  a) 


,  -ik;r(c0S    9  + cose)  , 

f^°(Tt,a,)f^°(9,n)  +  e       ^  f^°(o,a)f^^^(e,o) 


ik2d  r  ik-|(  cos  &  -  cos  a ) 


f^°(o,a)f^(n,o)f^°(9,Ti) 


-ik7?(cos9  -  cosa)  ,  , 

+  e  f     (n,a}f     (o,Ti;f     (9,o) 


ik3d 


ik'^(cos6+   cos  a)  ,  , 

2  ~bO/        \^ao/        \^bo/        \«ao/„      -> 

e  f     (n,a)f     (o,n)f     (n,o)f     (»,n) 


-ik--(  cos  9  +  cos  a ) 


ikd 


+  e 


ik^(cos9  +  cos  a) 


^30  /        % -bo  /        s .ao ,        V „bo  /  _      \ 
f     (o,a)f     (n,o)f     (o,n)f     (e,o) 


„bo,        x„2    -ao, 


<.bo, 


f^(TT,a)D^  f'^"(9,TT)+2D^f""(K,c)D^  f^°(9,:T) 


{jf^(n,c 


-ik7r(  cos  &  +  cos  a ) 


+e 


,a)+D|f^(n,a)Vf^(e,n) 
„a0/.     _  ^^2  ,bO/^      x.^   ^ao/        v„      „bo. 


^^  ^f^°(o,a)D|  f^(6,o)*2DQf^°(o,c)DQ  f^°(9,o)    \ 


Jf^(o,a)+D|f^°(o,a)|f^(e,o) 


^0'(d-^/2j 
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We  have  presented  in  (31)  a  relation  between  the  scattered  far  field 
of  the  combination  and  those  of  the  component  cylinders.  Vie   wish  to 
point  out  that  a  similar  relation  holds  between  the  corresponding  fields 
at  all  points  of  space.  But  this  relation  will  not  be  detailed  here. 

5.  A  special  case.  Scattering  by  two  conducting  circular  cylinders. 

The  abstract  relations  obtained  above  can  be  verified  in  the 
special  case  of  scattering  by  tvro  parallel  conducting  circular  cylinders 
A  and  B  with  radii  a  and  b  respectively.  If  the  electric  field  vector 
is  parallel  to  the  axis  of  the  cylinders,  we  are  confronted  with  the 

familiar  boundary  value  problem 

2 

D.S.   (^  +  k  )u  =  o 

(32)  B.C.  u  "  o  on  th^  boundaries  of  the  cylinders. 

(33)  R.C.         lijn    V^  (|2  -  ikU)   =  o 

r->oo 

Solution  of  this  problem  for  each  cylinder  in  isolation  requires  the 

determination  of  coefficients  a°  and  b  in  (3U)  and  (35)»  below 

n  n 

00 


(3U)  U°-i:     a°H;^^^(kr^)e^"^^ 


^n 
n=-oo 


and 


(35)  ^,=Z^:H<^>(Xry"<'^^' 


n=-oo 


Expansion  of  the  incident  plane  wave  in  terms  of  cylindrical  waves  by  means 
of  the  Fourier-Bessel 
conditions (32)  yields 


r9i 

of  the  Fourier-Bessel  expansion^-" and  application  of  the  boundary 
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00 


J  (ka)  /T>  in(&  -a) 

(36)     <f  ■  E      -i"  -JIT--    H<"(kr,)  e 
n=-co  H       (ka; 


n 


00  J  (kb)  /^N  in(0,-a) 

n 


for  the  unperturbed  scattered  fields,  as  is  well  known.     Using 
the  asymptotic  form  of  H         for  large  argument 

ikr       -ir 


H>.,~^e^/ir» 


we  see  that  the  far  fields  are 

^^^a       -iJI    /-5-  00       J  (ka)  in(a^^) 

(36)     K'^    e    ^/i       (-l)Z      -^r-T    ^         ' 
//r~  '  n=-oo  H'    '(ka) 

fa  n 


^^b       -ir    /2"  00        .  ,,,  X  in(e  -a) 

(39)     u°  =  e e    ^  / 15^       (-1)  ^I      V^^^         e        ^ 

)^  "=-°°  H^^^(kb) 

n 


Normalizing  our  coordinate  system  to  a  new  origin  midway  between 

the  origins  in  the  scatterers  as  we  did  above  in  the  abstract  treatment, 

we  obtain 

ik(iH-r  cos  B)  .,  d  .n      /»—  r—  -r  n     \ 

^y.^rf'.^-^ e-^"^"e-^E^    (-1)1        '-^      e^"<*-«>. 

n 

and 

ik(r-  I  cos  e)       ^^     a     -iJi    /r                 00       J  (kb)  .    ,^_„  x 

(Ul)     if.« e       2  ,     HV^       (.1)^         n    ^xn(e-a)   , 

^  (/?  ^  n-oo  H^^^(kb) 
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from  which  we  observe  that  the  unperturbed  patterns  are  respectively 
-ir  /fr  00   J^(ka) 


(as)  C(e,»)  =  e  ^  /I  (-i)E_  -b-T   = 


in(e-a) 


n=-oo  H^^(ka) 
n 

n 

The  perturbed  patterns  can  be  obtained  by  the  following  procedure; 

We  take  the  multiple  scattering  viewpoint,  as  we  did  in  the 
abstract  case.  Each  excitation  elicits  a  response.  The  responses  to 
the  initial  plane  wave  have  the  form  (3U)  and  (35).  These  responses, 
in  turn,  excite  the  second  cylinder,  eliciting  new  responses.  Unlike 
the  procedure  of  the  abstract  treatment,  however,  these  excitations  are 
not  expressed  in  terms  of  plane  waves,  nor  are  the  responses  expressed 
in  an  asymptotic  form. 

At  all  stages  intermediate  between  the  initial  plane  wave  excitation 
and  the  evaluation  of  the  far  field  response  of  the  combination,  the 
excitations  will  have  the  i^eneral  form 


(U.)   f   c  ^f   (kr)  e^"* 


and  the  responses  will  have  the  general  form 


q=-oo  ^  ^ 


The  proper  adjustments  must  be  made  in  (hU)  and  (U5)  for  angles 
of  incidence  and  the  fields  must  be  expressed  in  their  proper  coordinate 
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systems.  The  coefficients  C  and  d  are  determined  successively  for 
each  scattering  by  the  boundary  conditions  on  the  cylinders  which  are, 
in  this  case,  that  the  sum  of  the  excitation  and  its  response  is,  for 
each  cylinder,  zero  on  the  boundary  of  that  cylinder.  The  application 
of  the  boundary  conditions  requires  that  each  excitation  be  represented 

in  the  same  coordinate  system  as  its  response.  This  is  accomplished  by 

[91 
the  application  of  addition  theorems  for  cylindrical  waves'--'.  The 

process  is  continued  until  the  desired  degree  of  scattering  is  obtained. 

-3/2 

Application  of  this  method  yields,  up  to   terms  of  degree  d         ,   the 

following  perturbed  patterns: 


(Equations  on  pages  22  and  23) 


Throughout  this   calculation  we   follow  the  procedure  of  THerslnr'- -' J  but, 

-3/2 
unlike  that  author,  we  do  not  neglect  any  terms  of  order  d         ,  in  computing 

the  final  form  of  the  result  for  large  d. 
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Calculations  of  a  similar  nature  have  been  carried  out  by 
Twersky'--'   for  two  identical  cylinders  with  more  general  boundary  conditions. 
Comparison  of  (li6)  and   (hi)  with  formula   (6)  of  Twersky's  paper  show  that 

the  non-interaction  terms  in  our  result  agree  with  Twersky's  first  tern 

-1/2 
and  our  terms  of  degree  d  agree  with  Twersky's  second  term.     Our 

term  of  degree  d"     agrees  with  the  third  term  in  Twersky's  paper 

and  also  with  an  equivalent  corresponding  term  in  a  report  L^J  by  Twersky 

on  which  his  paper  was  based.     Our  terms  of  degree  d     '     do  not  agree 

with  Twersky's  fourth  term,  but  we  should  not  exi^ect  such  agreement. 

Twersky's  terms  are  the  asymptotic  representations  of  successive 

orders  of  scattering,   that  is  successive  bounces.     Our  successive  powers 

-1/2 

of  d         ,   however,  do   not  correspond  to  successive  bounces.     Our  term 

of  degree  d  -"     ,  for  example,  contains  contributions  fmm  the  second 
bounce  as  well  as  the  fourth  bounce. 

The  procedure  used  above  for  obtaining  (U6)  and  (Ij?)   is  long  and 
tedious.     These  results  need  not  be  obtained  by  that  procedure.     The  use  of 
the  abstract  formulas   (25)  and  (26)  simplify  the  calculation  considerably 
and  should  do  the  same  in  any  other  case  where  the  scattering  problems 
for  the  component  cylinders  are  separable.     All  we  need  do  is  substitute 
(U2)  and   (1x3)  into   (25)  and  (26)   respectively.     The  results  follow 
immediately  and  agree  with  (U6)  and   (hi)  respectively, 

6,     An  Explanation  of  the  Result 

The  expression  (31)   for  the  scattered  far  field  of  the  combination  of 
two  cylinders  appears,   at  first  glance,   to  be  complicated.     However,   it  is 
not  as  abstruse  as  it  seems.     A  closer  examination  of  these  expressions 
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reveals  the  significance  of  these  various  terms  and  factors.  We  note  first 
that  the  terms  are  grouped  in  increasing  orders  of  accuracy.  We  note  also 
that  they  are  grcuped  in  pairs.   The  first  members  of  each  pair  represent 
fields  ultimately  scattered  by  cylinder  ^   while  the  second  members  of 
the  pairs  repi-esent  fields  ultimately  scattered  by  cylinder  B.  The  first 

pair  corresponds  to  single  scattering  while  the  other  pairs  correspond 

ilkxcos  9 
to  multiple  scattering.  The  factors  e   ^     represent  the  phase  differences 

for  the  scatterers  relative  to  the  point  of  observation  while  the  factors 
j,d 
iik-KCos  a 

e         take  into  account  the  phase  of  the  incident  wave  at  the  center 
of  a  scatterers  when  the  incident  wave  has  zero  phase  at  the  origin. 
The  f 's  containing  a  &  dependence  are  scattering  patterns  whereas  the 
f 's  containing  specific  values  for  6  are  excitation  factors  accumulated 

in  the  multiple  scattering.  We  note  also  that  the  factors  of  the  form 

iknd 
e    where  n  =  0,  1,  2,  3  refer  to  the  increase  ikd  in  the  phase  of  a 

wave  in  going  from  one  scatterer  to  another  and  that  n  signifies  the 

number  of  bounces. 

As  an  illustration  of  the  above  remarics,  let  us  consider  some  of 

eikr  ik*(cos  &  -  cos  a) 
the  terms  in  more  detail.  The  term  ■    e  ^  f  (&,o) 


would  occur  in  the  case  of  no  interaction,  i.e.  in  the  limit  of  infinite 

spacing.  The  factor  exp  — x-(cos  &)  takes  account  of  the  fact  that  the 

origin  is  not  at  the  center  of  A.  The  positive  sign  preceding  cos  &  shows 

that  the  scattered  wave  came  from,  cylinder  A,  The  negative  sign  in  the 

.,d 

exponent  of  the  factor  e       ^  shows  that  cylinder  A  received  the  milial 


excitation.     It  gives  the  phase  of  the     incident  wave  at  A,     The  term 

^ikr    ik5(cos  9  *  cos  a) 

«  ^(n,a)>,r) 


/? 
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differs  from  the  first  term  in  the  follovd.ng  respects.     The  positive 
sign  in  the  factor  exp(+  — sr-cos  o)  shows  that  cylinder  B  received  the 
original  excitation  as  does  the  factor  f,        \.     The  phase  factor  e 
represents  the  increase  of  the  phase  of  the  wave  in  going  from  B  to  A, 

We  shall  now  explain  a  typical  term  of  the  last  square  bracket. 

The  first  term  in  the  last  square  bracket  of  (31)   represents  the 
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scattering  by  /,  of  a  term  of  order  d    initially  scattered  by  B. 

ikd 
The  phase  factor  e    takes  account  of  the  travel  of  this  wave  from 

B  to  A. 

The  positive  sign  preceding  cos  a  shows  that  cylinder  B  was 
excited  initially.  The  differentiation  of  the  scattering  amplitudes 
shows  the  effect  of  a  higher  order  excitation  of  A  by  the  field  initially 
scattered  by  B,  since  the  higher  order  response  of  B,  (which  acts  as 
an  excitation  for  A),  is  re presentable,  near  A,  as  a  derivative  of  a 
plane  wave  with  respect  to  angle  of  incidence. 

The  explanation  of  further  terms  in  the  final  result  proceeds 
on  the  same  lines  as  the  explanations  given  above.  We  omit  these 
explanations  for  the  sake  of  brevity. 

7,  The  total  scattering  cross-section  of  a  combination  of  two  identical 
circular  cylinders. 

We  consider  a  plane  wave  normally  incident  on  a  pair  of  identical 
circular  cylinders.  The  circumstances  are  illustrated  in  Figure  6. 
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1 


FIGURE  6 

The  computation  of  the  total  scattering  cross  section  6^  of  the 
t*o  identical  circular  cylinders  is  facilitated  by  the  use  of  (30)  in 
conjunction  with  the  following  well  known  theorem'- ■'*  *^"^> 


(U8)      <5-  -  -  ^  Re  F(|,^) 


where  F(ft,a)  is  the  scattering  amplitude  (30)  of  the  combination. 
The  phase  factors  are  simplified  for  the  values  &  ■  x  ,  a  »  •»,  The 
fact  that  the  cylinders  are  identical  enables  us  to  write  f °  -  f   «  f 

Further  simplifications  result  from  the  geometrical  sjTmnetry  of  the 
problem,  namely 


f°(o,n)  -  f°(T,,o),  f°(|,o)  -  f°(o,|)  -  f°(n^)  -  f°(|,n) 
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The  total  scattering  cross  section  in  teriris  of  the  spacing  is  then 


(U9)     or  -  -  g  Re 


l/?nIE  e 


ikd 


f°(c,^) 


2*  e 


ik2d 


f°(o,^) 


V(ti,o) 


ik3d 


f°(o,|) 


L°(Tt,0) 


ikd 
1   ^  ^  21k  (f°(o,|)D^°(o,|).2D,f°(o,|)D,^r5(o,|) 


L 


^f°(o,|).D2f°(o,|) 


+  0'(ci-^/2) 


f°(o,p 


8.     Additional  Remarks 
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It  is  clear  that  interactions  of  degree  greater  than  d  '" ''   can  be 
computed  by  the  inclusion  of  more  terins  in  the  expansions  used  to  obtain 
these  results. 

It  is  also  clear  that  the  method  applies  to  cases  of  more  than  two 
scatterers,  but  the  computations  would  be  more  tedious  than  in  the  case 
of  two  scatterers.  The  computations  might  be  simplified  by  using  a 
"  consistency"  method  employed  in  [2]  rather  than  tracing  the  successive 
scattering  in  detail.  This  method  involves  a  steady  state  point  of  view. 
The  response  of  each  cylinder  is  expanded  in  a  neighborhood  of  each  of 
the  other  cylinders.  Each  cylinder  will  then  be  excited  by  the  incident 
plane  wave  and  by  an  approximately  plane  wave  from  each  of  the  other  cylinders. 
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These  considerations  introduce  certain  undetermined  coefficients  which 
can  be  determined  by  imposing  the  requirement  that  the  field'  scattered 
by  the  various  cylinders  be  consistent  with  one  another.     EJvaluation  of 
the  coefficients  will  provide  the  solution. 
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II.  Three  Dimensional  Case. 


1,  Statement  of  the  Problem.  The  method  employed  previously  in  the 
case  of  multiple  scattering  of  plane  waves  by  two  widely  spaced 
cylinders  of  artitrary  shape  can  be  applied,  also  to  the  corresponding 
three  dimensional  scalar  problem  for  two  bodies  of  arbitrary  shape. 
The  assumptions  here,  as  in  the  previous  case,  are  that: 


1)  The  spacing  of  the  bodies  is  large  compared  to  their  dimensions. 

2)  The  spacing  is  large  compared  to  the  wavelength  of  the  incident 
plane  wave, 

3)  The  response  of  each  scatterer  to  a  plane  wave  excitation  is 
known  when  the  scatterer  is  isolated  in  space.  In  other  words,  the 
individual  unperturbed  complex  scattering  amplitudes  of  the  far  fields 
are  known. 

The  situation  is  the  following: 
A  plane  wave  of  unit  amplitude 

ik(x  sin  &  cos  ^^  +  7   sin  ©.  sin  ;6  +  z  cos  &  ) 
(1)  u  =  e 

where  &  and  (^    are  the  angles  of  incidence  (see  Figure  1),  is  incident 
upon  the  combination  of  two  bodies  (Figure  2,  page  32). 


Figure  1 
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The  response  of  each  body  in  isolation  to  the  incident  plane  wave  is  of 
the  form 


^ikr  00       f(^,»,^,(/) 

(2)       u  =v-i: 

3        r      ^ — 


n 

n=o  r 


Here  we  are  referring  to  spherical  coordinates,   r,  &,  'ji,  such  that 

(2a)  X  =  r  sin  9  cos  0 

y  =  r  sin  B  sin  flf 
z  =  r  cos  & 

The  letters  x,  y,  z  are  cartesian  coordinatesof  a  point  P  with 

respect  to  a  coordinate  system  located  appropriately.  The  quantities 

&  and  0  which  specify  the  direction  of  travel  of  the  incident  wave  can 
o     o 

be  thought  of  as  the  colatitude  and  longitude  respectively  of  a  point 

on  the  unit  sphere.  That  point  is  determined  as  the  intersection  with 

the  unit  sphere  of  a  line  passing  through  the  origin  in  the  direction 
of  travel  of  the  incident  wave. 

The  object  is  to  obtain  the  perturbed  complex  scattering  amplitude 
of  the  combination  as  a  function  of  the  individual  unperturbed  scattering 
amplitudes.     This   can  be  done,  as  in  the  previous  case,  by  first 
expanding  the  field  scattered  by  each  body  in  a  neighborhood  of  the 
other  body  and  then  expressing  these  scattered  fields  in  terms  of 
plane  waves  and  their  derivatives.     The  a  priori  assumption  that  the 
unperturbed  responses  to  an  incident  plane  wave  are  known  pennits  the 
calculation  of  the  perturbed  responses. 

We  must  establish  coordinate  systems  for  the  problem.     We  first 
circumscribe  spheres  A'  and  B'    about  the  respective  bodies  A  and  B, 
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This  establishes  the  origins  of  the  coordinate  systems  of  A  and  B 
(Figure  2) 


Incident 

Plane 
Wave 


FIGURE  2 


The  spacing  d  is  defined  as  the  line  joining  the  two  centers.  We  then 
choose  the  x  axis  (Figure  3).»»i«e.  the  x  axis  for  the  coordinate 
system  of  A,  to  be  coUinear  with  d  and  such  thct  x  increases  in  the 
direction  of  B» 


FIGURE  3 

The  y     and  z     axes  can  then  be  chosen  to  form  a  right-handed  system  and 
a  a 

the  coordinate  system  of  B  is  parallel  to  that  of  A.     (Choice  of  the 

z  axis  along  d,   instead  of  the  x  axis, would  be  more  symmetrical,  but  has 

been  avoided  for  reasons  which  will  be  explained  later.) 
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2.  Expansion  of  the  scattered  waves  in  terms  of  plane  waves » 

We  proceed  to  expand  the  response  of  each  body  in  a  neighborhood 
of  the  other  body.  This  exi^ansion  yields  the  following  field  scattered 
by  A  in  a  neighborhood  of  B  where  the  cartesian  coordinates  are  in  the 
B  system. 


(3)  U°  -  e"'('**"t; 


ao/n 


'^^C^v^o^'^V'  y^^C^i^^o^'^'V 


-3' 


<°(|,e,,o,0^)^*d(d-^) 


J 

In  calculating  equation  (3)  the  phase  of  the  incident  wave  at  the  'center' 
of  A  has  been  ignored.  The  appropriate  phase  factor  will  be  restored  in 
the  final  result.  A  similar  expansion  of  the  field  scattered  by  B  in  a 
neighborhood  of  A  yields 


(U)  U°  -  e^^^-^a^ 


d^o  %'^o'^'^o)*T 
d 


2.  2 


i!^^F^f^(i,«„.n,(^,)*vr(|.«o.",^o> 


-=^9^r<f»o""«'o)-5*0C(|'*o""l*o) 


*'l<--Z'%'"-V 


>■   + 


e'Cd"^) 


where  D^  -  §5  and  D^  =  ^ 


We  now  express  these  scattered  fields  in  terms  of  plane  waves  by 
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means  of  the  following  relations; 


-^^o'V 


=  e 


ik(x  sin  9     cos  0  +  y  sin  9     sin  0  +  z  cos  B   ) 
o  o  o  o  0    * 


v(i,0)         »  e 


ikx 


v^^(|,0) 


fH       \  -iKX 

,     v(^,n)  =  e  , 

-Ikze  ,    Vf.   ("5*")  "  -ikze 


o 


v^   (|,0)     =  ikye^      ,     vj^  q,n)  -  -ikye"^ 


1 


1 
1^ 


^9  e  (|'°)*  -fi  0  ^1>'^ 


o  o 


o  o 


0  0  oo 


^- 


2 


ikx 


-ilex 


Substitution  of  these  relations  into   (3)  and  (U)  yields: 


(5)    0^    -e 


Ikd 


KfiK^^>%fiA>  ^(^ 


"^sffo'-^^*!'"' 


o  o 


C(l'%.°''«o' 


o 


*T(|,0)fJ°(|,9^,0,<!f^)l.|3'(d-3) 


and 
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(6)  if  =  e 


Ucd 


|v(^,.)f^(^,9„,<H  ij  U^ 


Vq  q   (^,n)+V0  ^  (-jjrO 


o  o 


c^o 


C^i'%--V 


o 


These  fields  can  then  be  expressed  in  terms  of  f  alone  by  means  of  a  recursion 
formula  (7)  obtained  by  Sommerfeld  L-^J , 


(7)  2iic(n+l)fj^3^ 


"(-"^i^^f'i-^a^)*;^^ 


n 


(8) 


In  our  case  n  ■  0,  and  the  recursion  reduces  to 


2ikf  j^  -  .  cot  9  Dg  +  d|  +  csc^  ^0    f    ^o      • 


It  is  here  that  the  proper  choice  of  a  coordinate  system  is  seen 
to  be  convenient.    We  note  that  the  recvirsion  involves  the  Beltrami 
Operator  which  is  singular  at  9  -  0.     This  singularity  occurs  in  the 
direction  joining  A  to  B,  if  the  z  axis  is  chosen  collinear  with  d. 
Another  difficulty  in  using  such  a  coordinate  system  is  the  occurrence  of 
terms  in  the  expansions   (3)  and  (U)  which  are  not  easily  replaceable  by 
derivatives  of  plane  waves.     The  coordinate  system  we  have  selected 
overcomes  these  difficulties. 


For  &  "  7  and  0  =0,   (8)  becomes 
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(9)     fi(>o'°'^o)  =  ^\^l  -  4]  ^o^i>^o>^K^ 


For  ®  =  "s  and  ^  =  "*    (8)  becomes 


(ir))fl(|,^,,n,{^,)-2l 


2ik 


\*^0 


'    ^o(|'*2>"»^o^ 


Substitution  of  (9)  and  (IC)  into  (5)  and  (6)  results  in  expressions 
(11)  and   (12)  for  U°     and  U°   ,  which  do  not  contain  f^°and  f^    . 


(11)  U°     =  e 
s 
a 


ikd 


a 


i.(|.0)C(|.^.0.!'o)*  ^\ 


1 

2l5 


(d|.d|]v(|,o) 


(.ao/H 


C(i.»o'°.Co' 


(12)  0° 
s. 


ikd 


o 

*(y(r3)  I 


i.(|.w)f^(|,e„,n,?„)  *^  [^  [d|^.d|J  v(|,n) 

o 
♦  E^  v(|.n)D^f^(f  e„,.,!^„) 


C<l'%'"'^o' 


*  ^"'l-'^fH 


0(d-3) 


J 


[  f^(f  »„.n,0„) 
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3.  Calculation  of  the  perturbed  far  field  amplitude. 

The  method  employed  previously  in  the  two  dimensional  case  will 
now  yield,  up  to  interaction  terras  of  degree  d     ,   the  field  scattered 
by  the  two  bodies. 


The  final  result  mvist  take  into  account  the  difference  in  phase 
of  the  original  plane  wave  with  respect  to  the  two  scatterers  and  also 
the  difference  in  phase  of  the  scattered  fields  of  the  two  bodies  at 
the  point  of  observation.     This  can  be  done  by  a  normalization  of  the 
coordinate  systems  and  subsequent  consideration  of  the  phases  in  the 
normalized  coordinate  system.     By  a  normalized  coordinate  system  we  mean 
a  central  coordinate  system  whose  origin  lies  midway  between  the  origins 
for  A  and  B  and  whose  x,  y  and  z  axes  are  parallel  to  those  emanating 
from  the  origins  of  A  and  B,     If  we  let  a  and  p  denote  the  angles  of 
incidence  corresponding  to  colatitude  and  longitude  respectively  (Figure  U)* 
we  obtain  the  fbllowing  expression  for  the  field  scattered  by  the  two 
bodies*   (The  subscript   "o"  is  omitted  from  the  f's  for  convenience,) 
(Here  the  incident  wave  is  taken  to  have  zero  phase  at  the  origin  of  the 
normalized  system  of  coordinates.) 
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h.     An  Explanation  of  the  Result.  The  expression  (13)  for  the  scattered 
field  of  the  combination  of  two  bodies  appears,  at  first  glance,  to  be 
complicated.  However,  it  is  not  as  abstruse  as  it  seems.  A  closer 
examination  of  this  expression  along  the  lines  followed  in  the  two  dimensional 

case  reveals  the  significance  of  the  various  terms  and  factors.  We  observe 

-2 
first  that  there  are  eight  terms  up  to  the  order  of  accuracy  d 

calculated  here,  and  we  shall  refer  to  them  by  number  in  the  order  of 

their  appearance  above.  The  first  two  terms  represent  single  scattering 

while  the  others  represent  multiple  scattering. 

The  functions  f  represent  the  complex  scattering  amplitude  of  the 
respective  bodies  with  respect  to  an  origin  in  the  center  of  the  relevant 
body,  when  the  incident  wave  has  zero  phase  at  that  center.  The  odd 
numbered  terms  represent  the  fields  ultimately  scattered  by  body  A  whereas 

the  even  numbered  terms  represent  fields  ultimately  scattered  by  B. 

±ik|sin  e  cos  (jf 
The  factors  e  represent  the  phase  differences  for 

the  scatterers  due  to  the  use  of  the  normalized  coordinate  system,  while 

tiltesin  a  cos  p 
the  factors  e  ^  take  into  account  the  phase  of  the  incident  wave 

at  the  center  of  the  respective  scstterers.  The  f 's  containing  a  B   and  0 

dependence  are  scattering  patterns  while  the  f 's  with  specific  values  for 

B  and  (jf  are  excitation  factors  accumulated  in  the  multiple  scattering.  We 

note,  finally,  that  the  factors  of  the  form  e  where  n  ■  0,1,2  refer  to 

the  increase  ikd  in  the  phase  of  a  wave  in  going  fir>m  one  scatterer  to 

another  and  that  n  signifies  the  number  of  bounces.     If  one  uses  the  above 

comments  as  a  guide,  all  the  terras  can  be  explained  in  a  say  similar  to 

the  explanation,   given  earlier,  of  the  final  expression  for  the  scattered 

field  in  the  two  dimensional  case* 
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Attn:  Mr.  0.  C.  Woodyard 

U.  S.  Array  Signal  Engineering  Labs. 

Evans  Signal  Laboratory 

Beljnar,  New  Jersey 

Attn:  Technical  Document  Center 

Massachusetts  Institute  of  Technology 
Signal  Corps  Liaison  Officer 
Cambridge  39,  Mass. 
Attn:  A.  D.  Bedrosian,  Room  26-131 


Commander 
Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio  (10 
Attn:  N.  Draganjac,  WCLN(J-l4 

Conmandar 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,   Ohio       (2 
Attni     Mr.    Paul  Springer.   WCLRE-5 

Commander 

Air  Technical  Intelligence  Center 
Wright-Patterson  Air  Force  Base,   Ohio 
Attn:     AFCIN-tola 

Commanier 

Rome  Air  Development  Center 
Griffiss  Air  Force  Base,   New  York 
Attn:     RCSSTL-1 

Commander 

Rone  Air  Development  Center 
Griffiss  Air  Force  Base,  New  York 
Attn:  Mr.  Donald  Dakan,  RCOE 


Commanding  General,  SIOFM/EL-PC 
U.  S.  Armv  Signal  Engineering  Labs. 
Fort  Monmouth,  New  Jersey 
Attn;  Dr.  Horst  H.  Kedesdy 

Deputy  Chief,  Chem-Physics  Branch 

Commareier 

Army  Rocket  and  Guided  Missile  Agency 

Redstone  Arsenal,  Alabama 

Attn:  Technical  Library,  ORDXR-OTL 

Commanding  General 

U.  S.  Array  Signal  Engineering  Labs. 

Fort  Monmouth,  New  Jersey 

Attn:  SIGFM/EL-AT 

Department  of  the  Army 

Office  of  the  Chief  Signal  Officer 

Washington  25,  D.  C. 

Attn:  SIGRD-7 

Office  of  Chief  Signal  Officer 
Engineering  and  Technical  Division 
Washington  25,  D.  C. 
Attn:  SIGNET-5 

Guided  Missile  Fuze  Library 
Diamond  Ordnance  Fuze  Laboratories 
Washington  25,  D.  C. 
Attn:  R.  D.  Hatcher,  Chief  Microwave 
Development  Section 


Commander 
Rone  Air  Development  Center  (ARDC) 
Griffiss  Air  Force  Base,  New  York 
Attn:  Dr.  John  S.  Burgess,  RCE 

Commander 

Air  Force  Missile  Development  Center 
Holloman  Air  Force  Base,  New  Mexico 
Attn:  HDOIL,  Technical  Library 

Director  (2) 

n.  S,  Amy  Ordnance 

Ballistic  Research  Laboratories 

Aberdeen  Proving  Ground,  Maryland 

Attn:  Ballistic  Measurements  Laboratory 

Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  Maryland 
Attn:  Technical  Information  Branch 


)Arraed  Services  Technical  Information 

Agency 
Arlington  Hall  Station 
Arlington  12,  Virginia 

)  Library- 
Boulder  Laboratories 
National  Bureau  of  Standards 
Boulder,  Colorado 

National  Bureau  of  Standards 

Department  of  Commerce 
Washington  25,  D.  C, 
Attn:  Mr.  A.  G.  McNlsh 

National  Bureau  of  Standards 

Department  of  Commerce 

Washington  25,  D.  C. 

Attn;  Gustave  Shapiro,  Chief 

Engineering  Electronics  Section 
Electricity  and  Electronics  Div, 

Office  of  Technical  Services 
Department  of  Commerce 
Washington  25,  D.  C. 
Attn;  Technical  Reports  Section 
(Unclassified  only) 

Director 

National  Security  Agency 
Washington  25,  D.  C. 
Attn:  rA>  (331) 

Ho.  Air  Force  Cambridge  Research  Center 

Laurence  G.  Hanscom  Field 

Bedford,  Mass. 

Attn:  CROTLR-2     -  P.  Condon 


(5)Hq.  Air  Force  Cambridge  Research  Center 
Laurence  G.  Hanscom  Field 
Bedford,  Mass. 
Attn:  CROTLS  -  J.  Armstrong 


(2) 


(5)  Hq.  Air  Force  Cambridge  Research  Center 
Laurence  G.  Hanscom  Field 
Bedford,  Mass. 
Attn:  CRRD 


Director,  Avionics  Division  (AV) 
Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D.  C. 

Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washin^on  25,  D.  C. 
Attn:  Mr.  E.  Johnston,  Code  R33E 

Commander 

0.  S.   Naval  Air  Missile  Test  Center 

Point  Mugu,   California 

Attn:      Code  366 

U.  S.  Naval  Ordnance  Laboratory 
White  Oak 

Silver  Spring  19,  Maryland 
Attn:  The  Library 

Commanier 

U.  S.  Naval  Ordnance  Test  Station 

China  Lake,  California 

Attn;  Code  753 

Librarian 

U.  S.  Naval  Postgraduate  School 

Monterey,  Callfoi-nia 

Air  Force  Development  Field  Representative 
Naval  Research  Laboratory 
Washington  25,  D.  C. 
Attn:  Code  1072 

Director 

D.  S.  Naval  Research  Laboratory 

Washington  25,  D.  C. 

Attn:  Code  2027 

Dr.  J.  I.  Bohnert,  Code  5210 

D.  S.  Naval  Resesrch  Laboratory 

Washington  25,  D.  C.  (Dnclasslf led  only} 

Classified  to  be  sent  to: 

Director 

U.  S.  Naval  Hesearch  Laboratory 

Attn;  Code  5200 

Washington  25,  D.  C. 

Commanding  Officer  and  Director 

U.  S.  Nav>*  Underwater  Sound  Laboratory 

Fort  Trumbull,  New  London,  Cormecticut 

Chief  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D.  C. 
Attn:  Code  U27 

Commanding  Officer  and  Director 

U.  S.  Navy  Electronics  Laboratory  (Library) 

San  Diego  52,  California 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Kashingtc-n  25,  D.  C. 
Attn;  Code  Ad3 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Surface  Guided  Missile  Branch 
Washington  25,  D.  C. 
Attn:  Code  ReSl-e 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 
Attn:  Fire  Control  Branch  (ReSt) 

Department  of  the  Navy 

Bureau  of  Aeronautics 

Technical  Data  Division,  Code  lil06 

Washington  25,  D.  C, 

Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 
Attn:  Code  817B 


Commanding  Officer 

tl.  S.  Naval  Air  Developnent  Center 

Johnsvllle,  Pennsylvania 

Attn:  NADC  Library 

Commander 

r.  S.  Naval  Air  Teat  Center 
Patuxent  River,  Maryland 
Attn:  ET-31S,  Antenna  Branch 

Director 

Naval  Ordnance  Laboratory 

Corona,  California 

Cormanding  Officer 

T'.  S.  Naval  Ordnance  Laboratory 

Corona,  California 

Attn:  Mr.  W.  Horensteln,  Division  72 

Airborne  Instruments  Laboratory,  Inc. 

160  Old  Country  Road 

Klneola,  New  York 

Attn:  Dr.  E.  G.  Fubinl,  Director 

Research  and  Engineering  Division 

Alrcom,  Inc. 
}9i   Main  Street 
Wlnthrop,  Mass. 

American  Machine  and  Foundry  Company 

Electronics  Division 

lOR?  Commonwealth  Avenue 

Boston  1?,  Mass. 

^ttn:  Mrs.  Rita  Moravcsik,  Librarian 

Andrew  Alford,  Consulting  Engineers 
299  Atlantic  Avenue 
Boston  10,  Mass. 

Avion  Division 
ACF  Industries,  Inc. 
ROO  No.  Pitt  Street 
Alexandria,  Virginia 
Attn:  Library 

Battelle  Memorial  Institute 

?0?  King  Avpnue 

Attn:  Wayne  E.  Rife,  Project  Leader 

Electrical  Engineering  Division 

Columbua  1,  Ohio 

Bell  Aircraft  Corporation 

Post  Office  Pox  One 

Buffalo  5,  New  York 

Attn:  Eunice  P.  Hazelton,  Librarian 

Bell  Telephone  Laboratories,  Inc. 

Whlppany  Laboratory 

Whippany,  New  Jersey 

Attn:  Technical  Information  Library 

Pacific  Division 

Bendlx  Aviation  Corporation 

11^00  Sherman  Way 

North  Hollywood,  California 

Engineering  Library 

Attn:  Peggie  Robinson,  Librarian 

Bendlx  Radio  Division 

Bendlx  Aviation  Corp. 

E.  Joppa  Road 

Towson  li,  Maryland 

Attn:  Dr.  D.  M.  Allison,  Jr. 

Director  Engineering  and  Research 

Boeing  Airplane  Company 

Pilotless  Aircraft  Division 

P.O.  Box  3707 

Seattle  2(j,  Washington 

Attn:  R.P.  Barber,  Library  Supervisor 

Boeing  Airplane  Company 

Wichita  Division  Engineering  Library 

Wichita  1,  Kansas 

Attn:  Kenneth  C.  Knight,  Librarian 

Boeing  Airplane  Company 

Seattle  Division 

Seattle  llj,  Washington 

Attn:  E.T,  Allen,  Library-  Supervisor 

BJorksten  Research  Labs,  Inc. 

P.  0.  Box  26? 

Madison,  Wisconsin 

Attn:  Mrs.  Fern  B.  .Korsgard 


Convair,  A  Division  of  General  Dynamics 

Corp. 
Fort  Worth,  Texas 

Attn:   K.G.  Brown,  Division  Research 
Librarian 

Convair,  A  Division  of  General  Dynamics 

Corp. 
San  Diego  12,  California 
Attn:  Mrs.  Dora  B.  Burke, 

Engineering  Librarian 

Cornell  Aeronautical  Laboratory,  Inc. 
IJi?^  Genesee  Street 
Buffalo  21,  New  York 
Attn:  Librarian 

Dalmo  Victor  Company 
A  Division  of  Textron,  Inc. 
1515  Industrial  Way 
Belmont,  California 
Attn:  Mary  Ellen  Addems, 
Technical  Librarian 

Dome  and  Margolin,  Inc. 
29  New  York  Avenue 
Westbury,  Long  Island,  N.  Y. 

Douglas  Aircraft  Company,  Inc. 

P.O.  Box  200 

Long  Beach  1,  California 

Attn:  Engineering  Library  (C-250) 

Douglas  Aircraft  Co.,  Inc. 
f27  Lapham  Street 
El  Segundo,  California 
Attn:  Engineering  Library 

Douglas  Aircraft  Companv,  Inc. 
3000  Ocean  Park  Boulevard 
Santa  Monica,  California 
Attni   P.T.  Cline 

Eq.  Sec.  Reference  Files, 

Eq.  Eng.  A250 

Douglas  Aircraft  Company,  Inc. 
2"00  North  Memorial  Drive 
Tulsa,  Oklahoma 
Attn:  Engineering  Library,  D-2S0 

Electronics  Communication,  Inc." 
1830  York  Road 
Tlmonlnm,  Maryland 

Emerson  and  Cuming,  Inc. 
R69  Washington  Street 
Canton,  Mass. 
Attn:  Mr.  W.  Cuming 

Emerson  Electric  Kfg.  Co. 

8100  West  Florissant  Avenue 

St.  Louis  21,  Missouri 

Attn:  Mr.  E.R.  Breslln,  Librarian 

Sylvania  Elec.  Prod.  Inc. 
Electronic  Defense  Laboratory 
P.O.  Box  205  -  (Unci) 
Mountain  View,  California 
Attn:  Library 

Fairchlld  Aircraft  Division 
Fairchlid  Eng.  and  Airplane  Corp. 
Hagerstown,  Maryland 
Attn:  Library 

Farnsworth  Electronics  Company 
3700  East  Pontiac  Street 
Fort  Wayne  1,  Indiana 
Attn:  Technical  Library 

Federal  Telecommunication  Labs. 
?00  Washington  Avenue 
Nutley  10,  New  Jersey 
Attn:  Technical  Library 

The  Gabriel  Electronic* 
Division  of  the  Gabriel  Company 
135  Crescent  Road 
Needham  rieights  9U,  Mass. 
Attn:  Mr.  Steven  Galagan 


General  Electric  Advanced  Electronics  Center 
Cornell  University 
Ithaca,  New  York 
Attn:  J.  B.  Travis 

General  Electric  Company 

Electronics  Park 

Syracuse,  New  York 

Attn:  Documents  Library,  B.  Fletcher 

Building  3-l!i3A 

General  Precision  Laboratory,  Inc. 

63  Bedford  Road 

Pleasantvllle,  New  York 

Attn:  Mrs.  Mary  G.  Herbst,  Librarian 

OKiodyear  Aircraft  Corp. 

1210  Hassillon  Road 

Akron  15,  Ohio 

Attn:  Library  DA20  Plant  A 

Granger  Associates 

Electronic  Systems 

966  Commercial  Street 

Palo  Alto,  California 

Attn:  John  V.  N.  Granger,  President 

Grumman  Aircraft  Engineering  Corporation 
Bethpage,  Long  Island,  N.  Y. 
Attn:  Mrs.  A.  M.  Gray,  Librarian 

Engineering  Library,  Plant  No.  5 

The  Halllcrafters  Company 

UliOl  West  5th  Avenue 

Chicago  2(j,  Illinois 

Attn:  La Verne  LaGloia,  Librarian 

Hoffman  Laboratories,  Inc. 
3761  South  Hill  Street 
Los  Angeles  7,  California 
Attn:  Engineering  Library 

Hughes  Aircraft  Company 
Antenna  Department 
Microwave  Laboratory 
Building  12,  Room  2617 
Culver  City,  California 
Attn:  M.  b.  Adcock 

Hughes  Aircraft  Company 
Florence  and  Teale  Streets 
Culver  City,  California 

Attn:  Dr.  L.C.  Van  Atta,  Associate  Director 
Research  Labs. 

Hycon  Eastern,  Inc. 
7?  Cambridge  Parkway 
Cambridge,  Mass. 
Attn:  Mrs.  Lois  Seulowitz 
Technical  Librarian 

International  Business  Machines  Corp. 

Military  Products  Division 

590  Madison  Avenue 

New  York  33,  New  York 

Attn:  Mr.  C.F.  McElwaln,  General  Manager 

International  Business  Machines  Corp. 
Military  Products  Division 
Owego,  New  York 

Attn:  Mr.  D.  I.  Marr,  Librarian 
Department  U59 

International  Resistance  Company 
UOl  N.  Broad  Street 
Philadelphia  8,  Pa. 
Attn:  Research  Library 

Jansky  and  Bailey,  Inc. 
1339  Wisconsin  Avenue,  N.  W. 
Washington  7,  D.  C. 
Attn:  Mr.  Delmer  C.  Ports 

Dr.  Henry  Jaslk,  Consulting  Engineer 
298  Shames  Drive 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 

Electromagnetic  Research  Corporation 
711  ll4th  Street,  N.  W, 
Washington  5i  D.  C. 
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Lockheed  Aircraft  Corporation 

2555  N.  Hollywood  Way 

California  Division  Engineering  Library 

Department  72-75,  Plant  A-1,  Bldg.  63-1 

Bnrbank,  California 

Attnt  N.  C.  Harnois 

The  Martin  Company 

P.  0.  Box  179 

Denver  1,  Colorado 

Attn!   Mr.  Jack  McComick 

The  Glenn  L.  Martin  Company 
Baltimore  3,  Maryland 
Attnt  Engineering  Library 
Antenna  Design  Group 

Maryland  Electronic  Manufacturing  Corp. 
5009  Calvert  Road 
College  Park,  Maryland 
Attn:  Mr.  H.  Warren  Cooper 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,  Rhode  Island 

The  W.  L.  Maxson  Corporation 

U60  West  llith  Street 

New  York,  N.  Y. 

Attn:     Miss  Dorothy  Clark 

McDonnell  Aircraft  Corporation 

Lambert  Sarlnt-Louis  Municipal  Airport 

Box  516,    St.   Louis  3,   Missouri 

Attn!     R.  D.  Drtrich,   Engineering  Library 


McMillan  Laboratory,    Inc, 
Erownvllle  Avenue 
Ipswich,    Massachusetts 
Attn:      Security  Officer, 


Document  Room 


Melpar,  Inc.  (2) 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

Attn:  Engineering  Technical  Library 

Microwave  Development  Laboratory 
90  Broad  Street 
Babson  Park  57,  Massachusetts 
Attn:  N.  Tucker,  General  Manager 

Microwave  Radiation  Company  Inc. 

19223  South  Hamilton  Street 

Gardena,  California 

Attn:  Mr.  Morris  J.  Ehrlich,  President 

Chance  Vought  Aircraft,  Ino. 
93m  West  Jefferson  Street 
Dallas,  Texas 
Attn:   Mr.  H.  S.  White,  Librarian 

Northrop  Aircraft,  Inc. 
Hawthorne,  California 

Attn:  Mr.  E.  A.  Freltas,  Library  DeptJlU5 
1001  E.  Broadway 

Remington  Rand  tiniv.  -  Division  of  Sperry 

Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,  Pennsylvania 
Attn:  Mr.  John  F.  McCarthy 

R  and  D  Sales  and  Contracts 

North  American  Aviation,  Inc. 
1221ii  Lakewood  Boulevard 
Downey,  California 
Attn:  Engineering  Library  l»9S-115 

North  American  Aviation,  Inc. 
Los  Angeles  International  Airport 
Los  Angeles  1:5,  California 
Attn:   Engineering  Technical  File 

Page  Communirations  Engineers,  Inc. 
710  Fourteenth  Street,  Northwest 
Washington  5,  D.  C. 
Attn;   Librarian 

Phllco  Corporation  Research  Division 

Branch  Library 

Ii700  Wissachlckon  Avenue 

Philadelphia  hh,   Pa. 

Attn:  Mrs.  Dorothy  S.  Collins 


Pickard  and  Bums,  Inc, 

2lj0  Highland  Avenue 

Needhau  9U,  Mass. 

Attn:   Dr.  J.  T.  DeBettencourt 

Polytechnic  Research  and  Development 

Company,  Inc, 

202  Tillary  Street 

Brooklyn  1,  New  York 

Attn:  Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,  Mass. 
Attn:  Dr,  John  Ruze 

Radiation,  Inc. 

P.  0.  Drawer  37 

Melbourne,  Florida 

Attn:  Technical  Library,  Mr.  M.L.  Cox 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  York 

Attn:  P.  S.  Carter,  Lab.  Library 

HCA  Laboratories 
David  Sarnoff  Research  Center 
Princeton,  New  Jersey 
Attn:   Miss  Fern  Cloak,  Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electronic  Products 
Building  10,  Floor  7 
Camden  2,  New  Jersey 
Attn:   Mr.  Harold  J.  Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

Admi  nlstrator 

The  Ramo-Wooldridge  Corporation 
P.O.  Box  !t5lj53  Airport  Station 
Los  Angeles  U5,  California 
Attn:   Margaret  C.  Whitnah, 
Chief  Librarian 

Hoover  Micnjwave  Co, 
9592  Baltimore  Avenue 
College  Park,  Marvland 

Director,  USAF  Project  RAND 

Via:  Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Rantec  Corporation 

Calabasas,  California 

Attn:   Grace  Keener,  Office  Manager 

Raytheon  Manufacturing  Company 
Missile  Systems  Division 
Bedford,  Mass, 
Attn:   Mr,  Irving  Goldstein 

Raytheon  Manufacturing  Company 
Wayland  Laboratory,  State  Road 
Wayland,  Mass. 
Attn:   Mr.  Robert  Borts 

Raytheon  Manufacturing  Comrany 
Wayland  Laboratory 
Wayland,  Mass. 

Attn:  Miss  Alice  G.  Anderson, 
Librarian 

Republic  Aviation  Corporation 
Farmlngdale,  Long  Island,  N.  Y. 
Attn:   Engineering  Library 

Thru:   Air  Force  Plant  Representative 
Republic  Aviation  Corp. 
Farmlngdale,  Long  Island,  N.Y, 

Rheem  Manufacturing  Conpany 
9236  East  Hall  Road 
Downey,  California 
Attn:   J,  C.  Joerger 

Trans-Tech,  Inc. 
P.  0.  Hox   3I46 
Frederick,  Maryland 


Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  Calirornia 
Attn:   Librarv  -  unclassified 

Sage  Laboratories 
159  Linden  Street 
Wellesley  81,  Mass. 

Sanders  Associates 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.  R.  Wild,  Library 

Sandia  Corporation,  Sandla  3ase 

P.O.  Box  5''00,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroscope  Company 

Great  Neck,  Long  Island,  New  York 

Attn:   Florence  W.  Tumbull,  Engr.  Librarian 

Stanford  Research  Institute 

Menlo  Park,  California 

Attn:  Library,  Engineering  Division 

Sylvanla  Electric  Products,  Inc. 
100  s'irst  Avenue 
Waltham  5li,  Mass. 

Attn:  Charles  A.  Thornhill,  Report  Librarian 
Walthan  Laboratories  Library 

Systems  Laboratories  Corporation 
1I1P52  Ventura  Boulevard 
Sherman  Oaks,  California 
Attn:  Donald  L.  Margerum 

TRG,  Inc. 

17  Union  Square  West 

New  York  3,  N.  Y. 

Attn:  M.  L.  Henderson,  Librarian 

A.  S.  Thomas,  Inc. 

161  Devonshire  Street 

Boston  10,  Mass. 

Attn:  A,  S.  Thomas,  President 

Bell  Telephone  Laboratories 
Murray  Hill 
New  Jersey 

Chu  Associates 

".  0.  Box  387 
Whitcomb  Avenue 
Littleton,  Mass. 

Microwave  Associates,  Inc. 
Burlington,  Mass. 

Raytheon  Manufacturing  Companv 
Missile  Division 
Hartwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  Laboratory 
225  Crescent  Street 
Waltham,  Mass. 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  50U,  Sunnyvale,  California 
Attn:   Miss  Eva  Lou  Robertson, 
Chief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:   Dr.  W.  C.  Hoffman 

Commander 

AF  Office  of  Scientific  Research 

Air  Research  and  Development  Com.iand 

lUth  Street  and  Constitution  Avenue 

Washington,  D.  C, 

Attn:   Mr.  Ottlng,  SRY 

Westinghouse  Electric  Corp. 
Electronics  Division 
Friendship  Int'l  Airport  Box  7U6 
Baltimore  3,  Marvland 
Attn:  Engineering  Library 
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Wheeler  Laboratories,  Inc. 
122  Cutter  Mill  Road 
Oreat  Neck,  New  York 
Attn:  Mr.  Harold  A.  Wheeler 

Zenith  Plastics  Co. 
Box  91 

Oardena,  California 
Attn!  Mr.  S.  S.  Oleesky 

Library  Geophysical  Institute 
or  the  University  of  Alaska 
College 
Alaska 

University  of  California 

Berkeley  U,   California 

Attn:  Dr.  Samuel  Silver, 

Prof.  Engineering  Science 
Division  of  Elec.  Eng.  Electronics 
Research  Lab. 

University  of  California 
Electronics  Research  Lab, 
332  Corv  Hall 
Berkeley  h,   California 
Attn:  J.  R.  Whinnery 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
!j?00  Oak  'Srove  Drive 
Pasadena,  California 
Attn:  Mr.  I.  E.  Newlan 

California  Institute  of  Technology 
1201  E.  Callfomi  a  Street 
Pasadena,  California 
Attn:  Dr.  C.  Papas 

Carnegie  Institute  of  Technology, 
Sohenlev  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Prof.  A.  E.  Heins 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  G.  C.  Daljnan 

University  of  Florida 
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